AD-A277_537 


MARYLAND 

COLLEGE  PARK  CAMPUS 


STUDY  OF  SUPERCONVERGENCE  BY  A  COMPUTER-BASED  APPROACH. 

SUPERCONVERGENCE  OF  THE  GRADIENT  IN  FINITE  ELEMENT 
SOLUTIONS  OF  LAPLACE*S  AND  POISSON’S  EQUATIONS 


by 


L  BabuSka 
T.  Strouboulis 
C  S.  Upadhyay 
S.  K.  Gangar^ 


torn 

com 

CMS 

Technical  Note  BN-1155 

C7>B 

og 

and 

CMC  Report  No.  93-07 

Texas  Engineering  Experiment  Station 

The  Texas  A&M  University  System 

November  1993 


INSTITUTE  EOR  PWYSICAL  SCIENCE 
AND  TECHNOLOGY 

Dnc  QJTACT7  iriLIACTED  I 


94  3  24  022 


SCCumry  CLAlSinCATlOH  this 


REPORT  DOCUMENTATION  PAGE 


1.  fkgpomr  HUMtCM 


READ  INSTRUCTIONS 
BEFORE  COUPLETTNC  FORU 


>.  Mcci^iCMT's  Catalog  mumscii 


Technical  Note  BN-1155 


4.  riTLC  rantf  JuAlirl*)  *■  TY^C  OF  MCPONT  A  ^CNlOO  COVCMCO 

Study  of  Superconvergence  by  a  Computer-Based  Final  Life  of  Contract 
Approach . 

Superconv;rgence  of  the  Gradient  in  Finite  Elem(!y.C;»FONN.HO  ono.  refoat  ,.umokn 
Solutions  of  Laplace  s  and  Poisson's  Equations 


7.  AuTMOAf.;  ^ 

1.  Babuskal-  T.  Strouboul is^-  C.  S.  Upadhyay^- 

S.  K.  Gangaraj2 


>■  fcnfonming  oncanization  HAMC  AMO  Aooness 


Institute  for  Physical  Science  and  Technology 
University  of  Maryland 
College  Park,  MD  20742-2431 


la.  mkooaam  eueucNr.  FMOjccT,  TASK 

AAKA  t  «OMK  UNIT  NUMOCHS 


IZ.  nC^ORT  OATl 


II.  controlling  orricc  namc  ano  aoorcss 

Department  of  the  Navy 
Office  of  Naval  Research 
Arlington,  VA  22217 


4.  MOniTORtNG  ACCnCY  NAMC  A  AOOnZtyU  4iiUfmt  C«icrW/iAA  Orttf)  if.  SC^UMITY  CL>SS.  (mi  iM«  fmpmn) 


November  1993 


IJ.  NUUSCR  OF  races 

59  +  Figures 


is«.  declassification/ OORNCRAOINC 
schedule 


IS.  OISTRISUTION  statement  (a/ tfiit  RApATtj 

Approved  for  public  release:  distribution  unlimited 


n.  distribution  statement  (•!  l/i«  aSacrao  aniararfin  RIacA  30,  Jf  ttUmrmt  fraai  XaRart; 


1*.  KEY  ROROS  rCani/nua  an  rnaaraa  a/Oa  It  nacaaaarr  anO  JOanM/F  Sr  »<*«*  niaa»ar> 


1 20.  abstract  r 


This  paper  addresses  the  problem  of  existence  of  the  superconvergence  points 
by  a  computer  based  proof.  We  prove  a  basic  mathematical  theorem  that  the 
superconvergence  point  exists  if  and  only  if  it  can  be  found  by  certain  numerical 
algorithm.  We  address  the  problem  of  the  superconvergence  points  for  the  gradient 
of  the  finite  element  solution  of  the  Laplace  and  Poisson  equations.  Our  study 
shows  that  the  sets  of  superconvergence  points  are  very  different  for  these  two  cases. 
We  consider  triangular  as  well  as  quadrilateral  elements  of  degree  p,  1  <  p  <  7. 
In  the  case  of  quadrilateral  elements  we  analyze,  among  others,  the  tensor-product 
and  the  serendipity  elements. 


coition  OF  I  HOV  ts  IS  obsolcts 
S  N  Oio:-  LF-  0'4-  4601 


security  cl ASSIFiCATion  of  TmiS  Rage  i’BA-.  Dmim 


study  of  Superconvergence  by  a  Computer-Based  Approach. 
Superconvergence  of  the  Gradient  in  Finite  Element 
Solutions  of  Laplace’s  and  Poisson’s  Equations 


I.  Babuska* 

Institute  for  Physical  Science  and  Technology  and  Department  of  Mathematics, 
University  of  Maryland,  College  Park,  MD  20742,  U.S.A. 


T.  Strouboulis^,  C.S.  Upadhyay^  and  S.K.  Gangaraj^ 

Department  of  Aerospace  Engineering,  Texas  A&M  University, 
CoUege  Station,  TX  77843,  U.S.A. 


November  1993 


*The  work  of  this  author  was  supported  by  the  U.S.  Office  of  Naval  Research  under  Contract 
N00014-90-J-1030  and  by  the  National  Science  Foundation  under  Grant  CCR-88-20279. 

iThe  work  of  these  authors  was  supported  by  the  U.S.  Army  Research  Office  under  Grant 
DAAL03-G-028,  by  the  National  Science  Foundation  under  Grant  MSS-9025110  and  by  the  Texas 
Advanced  Research  Program  under  Grant  TARP*71071. 


1 


Abstract: 

This  paper  addresses  the  problem  of  existence  of  the  superconvergence  points 
by  a  computer  based  proof.  We  prove  a  basic  mathematical  theorem  that  the 
superconvergence  point  exists  if  and  only  if  it  can  be  found  by  certain  numerical 
algorithm.  We  address  the  problem  of  the  superconvergence  points  for  the  gradient 
of  the  finite  element  solution  of  the  Laplace  and  Poisson  equations.  Our  study 
shows  that  the  sets  of  superconvergence  points  are  very  different  for  these  two  cases. 
We  consider  triangular  as  well  as  quadrilateral  elements  of  degree  p,  1  <  p  <  7. 
In  the  case  of  quadrilateral  elements  we  analyze,  among  others,  the  tensor-product 
and  the  serendipity  elements. 
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1  Introduction. 

Let  be  a  one  parameter  sequence  of  finite  element  solutions  of  a  problem 
which  are  computed  using  a  sequence  of  meshes  T  =  {T^}  and  let  u  denote  the 
exact  solution.  Let  us  assume  that  we  are  interested  in  the  values  of  the  solution 
or  its  derivatives  or  combinations  of  these  quantities  i.e.  in  the  linear  functional 
F(u)(as).  Let  us  assume  that  for  every  element  r  of  the  mesh  Tf^  a  special  point  0, 
which  depends  on  the  geometry  (but  not  the  size)  of  the  element,  is  given.  Then 
denoting 


(1.1)  $(tt  -  ttfc)  =  max  |F(tt  -  ttfc)(*)| 

we  are  interested  in  the  values  of 


(1.2) 


0(0;F;u,Ufc,h,T)  = 


\F{u-u^)ii)\ 


If  the  point  x  is  such  that 


(1.3)  Q{x\F;u,Uf^yhyr) — »0  as  h — >0 

then  X  will  be  called  a  u- superconvergence  point  relatively  to  the  exact  solution  u 
and  the  family  of  meshes  T.  Consider  now  a  family  U  of  solutions;  the  point  x  will 
be  a  U- superconvergence  point  if  it  is  «-superconvergence  point  for  every  u  €  U. 
Obviously  ©(i;  F; it, tt^, h, r )  <  1.  If  $(tt  —  u^,)  =  0  in  t  then  we  define 
0(sB;F;ii,it,„/i,T)  =  0.  Practically  we  can  be  interested  only  in  the  case  that 


(1.4)  0(*;F;ti,tt;„h,r)  <  q  <  1 

uniformly  for  all  it  6  2/  and  €  T.  In  the  present  paper  we  will  address  the  case  of 
superconvergence  in  the  sense  (1.3).  The  superconvergence  of  this  type  was  studied 
by  many  investigators;  see  for  example  [1-27]  and  the  dtations  in  these  papers.  The 
17-superconvergence  was  not  analyzed  in  the  literature;  we  will  address  this  type 
of  superconvergence  for  meshes  of  triangles  (resp.  meshes  of  quadrilaterals)  in  [28] 
(resp.  [29]). 

In  this  paper  we  will  study  the  superconvergence  problem  when  F(u)  is  a  deriva¬ 
tive  of  It.  Here  we  will  address  in  detail  only  the  cases  of  the  Laplace  and  the  Poisson 
equation  for  meshes  of  triangidar  and  quadrilateral  elements  of  degree  p.  We  will 
address  the  elasticity  problem  in  [30]. 
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Our  approach  is  completely  different  than  all  other  studies  of  superconvergence. 
In  particular  our  approach  is  computer-based  proof:  We  prove  a  theorem  that  the 
auperconvergence  point  exists  if  and  only  if  it  can  be  found  by  a  numerical  algorithm. 
Based  on  this  theorem  we  wrote  a  computer  program,  by  employing  this  program 
we  can  decide  whether  the  li-  or  u-superconvergence  points  exist  or  not  and  if  so 
where  they  are  located.  By  this  approach  we  can  handle  automatically  elements  of 
arbitrary  degree  and  2/-superconvergence  for  any  class  of  solutions.  This  is  essential 
because  often  in  practice  we  treat  only  homogeneous  differential  equations. 

Given  a  family  of  finite-element  meshes  and  a  corresponding  family  of  conform¬ 
ing  piecewise  polynomial  spaces  of  degree  p,  the  superconvergence  points  (for  a 
quantity  of  interest)  in  an  element  in  the  interior  of  the  domain  may  or  may  not 
exist  depending  on 

1.  The  geometry  of  the  grid. 

2.  The  class  of  solutions  of  interest. 

S.  The  span  of  the  local  polynomial  spaces  in  each  element. 

Following  this  Introduction  we  outline  the  model  problem  of  orthotropic  heat- 
conduction;  here  we  study  in  detail  the  superconvergence  for  the  isotropic  case 
(i.e.  Laplace’s  and  Poisson’s  equations)  but  the  same  approach  can  be  employed 
for  the  general  orthotropic  case.  We  present  the  theoretical  setting  of  the  study 
and  describe  the  computer-based  approach  which  gives  the  asymptotic  properties 
of  the  gradient  of  the  error  in  the  interior  of  the  grid.  We  then  demonstrate  that  the 
numerical  approach  can  be  employed  to  find  the  classical  superconvergence  points 
(if  they  exist)  in  the  interior  of  locally  periodic  finite  element  meshes  of  triangular 
or  square  elements  of  any  polynomial  degree. 

2  Preliminaries. 

2.1  The  model  problem 

Let  us  consider  the  model  problem  of  heat-conduction  in  orthotropic  medium 
with  mixed  boundary  conditions.  Let  Q  C  R’  be  a  bounded  polygonal  domain 
and  let  its  boundary  dtl  be  split  into  two  parts  Fd  and  Tff.  Let  u  be  the  solution 
of  the  problem 


(2.1o) 

(2.16) 

U  =  0 

on  Fd  , 
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(2.1c) 


3 

S  9*(«K  =  9  on  Ts  . 

h=l 

Here  Kuy  ft,  /  =  1,  2  axe  the  entries  of  the  ihermal-conduetivUy  matrix  which 
satisfy  the  conditions 


(2.2a) 


Kiu  =  Ktk  ,  kyl  =1,2 


(2.2i)  0  <  +  fl)  <  Y,  Klu(u(,<K^(il  +  (l)  Vt  =  ((„6)€R*, 

where  K^ux  ue  the  principal  thermal  conductivities; 

’  /Ul 

(2.3)  9*(«):=E^*r— .  ^  =  1.2 

are  the  components  of  the  flux  (heat-flux);  n*,  ft  =  1,  2  are  the  components  of  the 
unit  outer  vector  to  dH;  f  €  X^(fl),  §  €  L^Tif)  are  given  data. 

Remark  2.1.  The  majority  of  the  numerical  results  will  be  presented  for  the  case 
of  Poisson’s  equation  for  which  =  ,  where  6^  is  Kronecker’s  delta. 

Let  JIp^(f2)  :=  |ti  €  fr^(f2)  |  u  =  0  on  ri)|.  Then  the  above  problem  may  be 
put  in  the  variational  form:  Find  u  €  H^^(Q)  such  that 


(2.4)  Bn(uyv)  =  Ln(v)  V  u  e  ^^^(0) 

where 


(2.5a) 

and 

(2.56) 


•=  L  9v . 

Jn  JTif 
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If  Fp  =  0  we  assume  that  the  data  f^g  satisfy  the  consistency  eondUion 


(2.6) 


=  0 


We  define  the  enetyy-norm  by 


(2.7)  Il|«|lli  :=  B„(«.») 

Let  T  :=  {Th}  be  a  rtgular-famUy  of  meshes  of  triangles  or  quadrilaterals. 
(For  the  meshes  of  triangles  it  is  assumed  that  for  any  elements  r^,  Tj  6  the 
intosection  rjDrj  is  other  empty,  a  votex  or  a  common  edge,  and  that  the  minimal 
angle  of  all  the  elements  is  bounded  below  by  a  positive  constant,  the  same  for  all 
meshes.  For  the  meshes  of  quadrilaterals  we  employ  similar  assumptions  as  in 
[8],  [9].  The  meshes  Tk  are  not  assumed  to  be  quasiuniform.)  We  introduce  the 
conforming  finite-element  spaces 

(2.8)  C"(fl) I  »!..  o €  S'(f ) ,  k  =  l,..., M(T»)} , 

where  Fr^  is  a  vector-valued  mapping  function  for  the  A;-th  finite  element  (for  the 
proper  assumptions  on  the  smoothness  of  the  transformation  see  [8],  [9],  [20],  [26]); 
^(f)  denotes  the  element-space  of  polynomials  defined  over  the  master-element 
f;  M{Th)  is  the  number  of  elements  in  the  mesh  Th. 

Below  we  will  consider  the  following  choices  for  the  element-space  ^(t): 

a.  Complete  polynomial  space  up  to  degree  p. 

In  the  case  of  triangular  elements  we  have  5p(f)  =  Vpir)  where 

(2.9)  %{f)  ;=  { P  I  P(x, .  X,)  =  £  xj  x5  } . 


In  the  case  of  quadrilaterals  we  will  consider  the  following  choices  for  the  defi¬ 
nition  of  the  polynomial  space  5p(f)  (see  also  [31]). 

6.  Tensor-product  (hi-p)  polynomial  space  of  degree  p. 

(2.10)  S<''Kf):={p\  P(x. ,  X,)  =  X;  Oij  xj  xi  } 

'  *.2  ■' 
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c.  Serendipity  (trvnc)  polynomial  space  of  degree  p. 


(2.11)  ^(f):=|p|  =  +  + 

d.  Intermediate  polynomial  space  of  degree  p. 


(2.12) 


'■  io.  h=0  * 


•<<+*<» 


Remark  2.S.  Let  f  denote  the  master  square-element.  We  note  that: 


(2.13a) 

5('-‘)(f)  =  5'(f)  =  5'*(f), 

(2.136) 

5’(f)  =  5'*(f ) 

while  for  p  >  3  all  three  spaces  are  distinct. 

We  let  To  ^  finite  element  solution  of  the  model 

problem  satisfies:  Find  €  ^,ro 


(2-14)  BQ{Uf^,v^^  )  —  ^  €  ‘^.Tc 

We  let  Cf^  :=  u  —  denote  the  error  of  this  approximation.  We  assume  that  one 
of  the  above  mentioned  spaces  is  employed  with  p  >  1. 

2.2  The  class  of  locally  periodic  meshes 

We  now  present  the  definition  of  a  special  class  of  locally  periodic  meshes.  In 
the  following  Section  we  will  prove  several  results  for  the  asymptotic  behavior  of 
the  error  for  these  grids.  Let  us  consider  a  locally  periodic  grid  defined  as  follows: 

Let  0  <  fT  <  ,  X®  =  (x®,x®)  €  fl , 

(2.15)  5(«®,B)  :=  I®  =  (xi,xj)|  |x<-x?|<B,  t=l,2| 
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and  assume  is  sufficiently  small  such  that  5(9**,  W*)  C  H.  Further,  let  7  be  a 
set  of  multi-indices  (t,  j) ,  =  (*1***^  *a  **^)  €  fl  and 

(2.16)  c(.W>.A);=5(.<«>,A)cS(«",Jr),  (i.j)  €  7 

be  the  set  of  the  h-cells  (or  cells)  which  cover  exactly  5(8°,  H)  i.e. 

(2.17a)  U  c(.W),A)=S(e“,ff) 

(2.176)  c(8(*“»>,h)nc(8(^-*>,h)  =  0  for  (»i,ii)  (ta.ja) 

We  will  refer  to  5(8°,  H)  as  the  subdomain  of  periodicity  of  the  mesh  centered  at 
8°.  Denoting  by 

(2.18)  c  :=  |(8x,®a)  |  l*i|  <  1 ,  l*al  <  l| 

the  umt-  (master-)  ceU  c,  the  h-cell  is  an  h-scaled  and  translated  master-cell. 

Let  T  be  a  triangular  mesh  on  the  master-cdl  (the  master-mesh)  and  be 
the  mesh  on  c(8^'>'l,h)  which  is  the  scaled  and  translated  image  of  T.  We  will 
consider  the  family  T  of  locally  periodic  meshes.  Let  Tf^^T  and  be  the 

restriction  of  on  5(8°,  if)  and  the  restriction  of  T^(m^,S)  on  c(m^*'*\h). 
We  assume  that  (»,j)  €  7  i.e.  is  made  by  the  periodic 

repetition  of  the  h-scaled  master  mesh.  An  example  of  such  a  mesh  is  shown  in 
Fig.  la  where  the  master-mesh  is  shown  in  Fig.  2(a).  Note  that  the  parameter  2h 
has  the  meaning  of  the  length  of  the  cell  side  in  5(8°,  H).  In  D  —  5(8°,  E)  it  has, 
in  general,  nothing  common  with  the  size  of  the  elements  used. 

Below  we  will  study  the  asymptotic  properties  of  the  error  in  the  locally  periodic 
mesh  for  h-cells  well  within  the  interior  of  5(8°,  J7).  By  asymptotic  properties  we 
mean  the  properties  in  the  limit  as  h  tends  to  zero.  Moreover,  it  will  be  assumed 
that  H  >  C7h° ,  with  0  <  a  <  1  i.e.  the  subdomain  of  periodicity  of  the  mesh  in 
the  neighborhood  of  8°,  5(8°,  H),  can  shrink  to  zero  but  at  a  slower  rate  thu  the 
size  h  of  the  cells.  Nevertheless,  for  simplicity,  we  will  assume  that  H  is  fixed  and 
independent  of  h.  We  underline  that  changing  h  does  not  change  the  geometry  of 
the  mesh  in  5(8°,  Jf). 


3  Theoretical  setting 

3.1  Preliminaries 

We  will  here  for  nmplidty  address  only  the  case  of  Poisson’s  equation  (Kj^  =  Ski 
in  (2.1)).  Nevertheless,  the  results  could  be  easily  generalised  to  the  general  setting 
for  orthotropic  medium. 

The  solution  of  the  problem  is  understood  in  the  sense  of  (2.14).  In  (2.7)  we 
defined  the  energy-norm  which,  for  K  unit-matrix,  is 


(3.1)  llltillln  =  II  \Vu\  lln 

where  ||  •  ||  is  the  usual  X’-norm.  The  fact  that  we  deal  in  some  cases  with  the 
seminorm  instead  of  the  norm  will  not  play  any  role.  For  a  specified  domain  5  C  H 
we  will  also  use  the  norms  (seminorms)  |||  *  |||5,  ||  •  Us  with  obvious  meaning  as 
well  as  the  notation  B5(tt,v)  with  |||u||||  :=  B5(«,«).  Given  a  function  u  and  the 
multi-index  a  :=  (ai,a))  we  define 


(3.2a) 


zr«  := 


dxf'dxf*  ’ 


|a|  :=  aj  -|-  Oj 


(3.26)  (^^«)(»)  :=  l^«l’)(»)]*  >  integer 

and  hence  we  can  write 


(3.3)  |||«|||J  =  E  ||C-«|||  =  II  |V«|  ||>  =  ||®'.||J  . 

M.I 

Further  we  denote  by  |  *  |  the  usual  X”°-norm  and  ||  u  ^  =  max  |  D^u  | . 

|o|=l 

3.2  The  basic  properties  of  the  finite  element  solution 

Let  XT’  >  >  Hi  >  Hi  >  XTs;  we  haven  D  5(»®,Xr®),  S{m^,H^)  D  S{m^,H), 

and  S  {m°,Hi)  D  S  (as®,  XT,,,)  ,l<m.  We  will  deal  with  5  («®,Xr£),/  =  1,2,3,  with 
various  XT/  <  XT.  We  will  assume  that  Hi  are  such  that 
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(3.4tt)  Sim\Ht)=  \J  ,  /=  1,2,3. 

where 

{3.4»)  7,  :=  C  . 

This  means  that  S  (a:*’.  Hi)  can  be  partitioned  into  the  set  of  cells.  In  Fig.  lb  we 
give  an  example  of  the  subdomains  5(0**,  J7/),  f  =  1,2,  for  the  mesh 

shown  in  Fig.  la. 

We  will  make  the  following  assumptions  about  the  exact  solution  u: 

Assumption  I 
On  5(»®,F) 

(3.5)  1  1  <  if  <  oo,  0  <  |o|  <  p  +  2 

Remark  S.l.  Assumption  I  states  that  the  solution  is  smooth  in  the  subdomain 
5(e°,fr)  i.e.  the  subdomain  should  be  suiHdently  far  from  boimdaries,  material- 
interfaces  and  points  where  the  data  are  rough. 

Assumption  II 

If  a„  ;=  (i?“«)  («"),  a  =  (Qi,aj),  |a|  <  p -|- 1  then 


(3.6)  R}=  53  a»>0 

l«l=p+* 

Further,  we  assume  that  the  mesh  T(0,h)  is  such  that: 

Assumption  III 

On  S  (*®,  J7j),  H^<H^<H<H^ 

(3-7)  l|cfclls(ir«.jr,)  <  Ch^Hi 

with  >  (p  +  1)  —  e,  where  0  <  e  <  1  will  be  specified  later,  and  where  C  is 
independent  of  T(0,h),  H2,  but  it  depends  on  K  and  R. 
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Remark  S.2.  We  do  not  assume  that  it  is  smooth  in  Q  outside  of  For 

example,  can  have  a  boundary  with  corners  (as  in  Fig.  1(a))  and  hence  u  can 
be  unsmooth  in  the  neighborhood  of  these  corners.  Nevertheless  assumption  III 
makes  an  implicit  requirement  on  the  (refinement  of  the)  mesh  in  the  neighborhood 
of  these  comers.  If  u  is  smooth  in  a  convex  f2  and  the  mesh  is  quasi-uniform  then 


(3.8)  I  Cfc  I  n  <  I  I  o  .  »*  >  0 

and  hence  in  the  assumption  III  we  can  take  fi  <  p+1  arbitrary. 

Remark  S.S.  Assumptions  II,  III  imply  that  the  principal  part  of  the  error  in 
5(c°,  H2)  is  related  to  the  non-zero  (p  -|-  l)-derivatives  of  the  exact  solution  at  s*’. 

In  the  following  we  prove  a  series  of  lemmas  and  the  theorems.  We  w^l  assume 
everywhere  below  that  is  sufficiently  small  {depending  on  K,  R)  independent 
of  h.  For  simplicity  we  will  address  explicitly  only  the  case  of  triangular  meshes 
although  the  arguments  are  general. 

Lemma  1. 


(3.9a)  ll|ej|l.(.(«l.»)  £  V  c(»W),4)  C 


(3.96)  II  >  C6'  Vc(.(«),6)CS(«“,H,) 


(3.10a)  lll®Allls(.«.Hj)  >  »— 1,2,3 


(3.106)  II  aj,  lls(..,ff,)  ^  *  t  —  1,2,3 

where  C  depends  only  on  Hi,  R,  K  and  the  master-cell  mesh,  but  is  independent 
of  h. 

Proof.  Using  assumptions  I  and  II  the  lemma  follows  immediately  from  the  fact 
that  TJ^^u  is  bounded  below  on  S{x°,Hi),  and  therefore  on  c(®(*-^),h),  by  a 
positive  constant,  which  depends  on  R,  Hi  and  K.  Q 

Remark  3.4-  The  lower  estimates  in  lemma  1,  as  given  above,  hold  only  for  the  case 
of  triangular  elements  of  degree  p  (the  span  of  the  shape-functions  in  each  element 
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includes  all  monomials  up  to  degree  p).  In  the  case  of  quadrilateral  elements  the 
span  of  the  element  shape-functions  is  larger  and  includes  monomials  of  degree 
greater  than  p  (see  [31]).  Hence  here  we  have  to  understand  lemma  1  in  the  sense 
that  (3.9a),  (3.9b),  (3.10a),  (3.10b)  are  valid  for  the  solutions  which  are  not  in  the 
span  of  the  element-space. 

Remark  3.5.  We  also  remark  that  in  (3.9b)  we  have 


(3.11) 
and  not 


max 


1^1 

\  *  dXi  *  ’ 

1  1 

de^ 


(3.12) 

Nevertheless  later  we  will  assume  that 


>Ch\  t  =  l,2. 

5aJ,-  I  e(B(*J),A) 


(3.13) 


'a®,. 


with  c  uniform  for  the  entire  class  of  functions  under  consideration.  This  can  be 
achieved  by  imposing  additional  assumptions  on  in  assumption  II. 

Lemma  2.  There  is  a  polynomial  Q  of  degree  p  1  on  5(®®,  ffj),  J?2  <  Hi 


(3.14a)  |||u  —  Qlll5(a0_jy,)  <  C 


(3.146)  I  ^“(1*  -  <3)  1  <  CAT’''”' ,  0<H<p  +  2 

Proof.  Q  is  the  Taylor  expansion  of  u  up  to  degree  p  -f  1.  □ 

Lemma  S.  Let  Q  be  a  polynomial  of  degree  p-f- 1  defined  on  S  (as®,  H)  and  let  Q™"* 
be  its  interpolant  of  degree  p  on  the  mesh  Th  (as®,  H).  Denote 

(3.15)  p:=Q-  Qr 
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then  p  is  a  periodic  function  on  S  (e°,  J7)  with  poriod  2h,  p  €  and 

hence 


(3.16) 


p(xS**’^  +  K *a)  =  -  K *3).  1*2  -  *3 <  h 

p(«i, +  k)  =  p(xi,xS,*‘'^  -  h),  |xi  -  <  h 


Proof.  The  lemma  follows  immediately  from  the  observation  that  the  (p  +  1)- 
derivatives  of  Q  are  constant  in  5  Note  that  we  are  dealing  with  elements 

which  are  able  to  reproduce  any  polynomial  of  degree  p  on  Q 

In  the  following  we  will  consider  polynomials  of  degree  (p  +  1) 


(3.17)  Q(xi,xa):=  6a  (®i  -  a5?)“‘(x,  -  x5)“*  ,  a  =  (ai,aa) 


and  we  assume  that 

(3.18)  0  <  pj  <  5;  6’  <  <  00 

|a|s|H-l 


Lemma  4.  Let  p  be  given  by  (3.15);  then 


(3.19a)  ,  J  =  0,  I 

(3.19»)  C, <\Vp\  ,  j  =  0,  1 

where  Ci,  Ca  depend  on  pi  and  pa  &nd  the  master  mesh  but  are  independent  of  h. 

Proof.  The  lemma  follows  from  the  basic  properties  of  the  interpolant,  (3.18)  and 
lemma  1.  Q 

We  now  define  the  periodic  finite  element  problem  on  c(«^**^\6).  Let 
(3.20)  :=  ju  G  :u  satisfies  (3.16)| 

and 
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(3.21)  (c(.W>,  h))  ;=  €  Bl„  (<<.«••■>,  h))  ■.  u  |,€  t  €  li«'| 

Let  ^  €  Snm  denote  z*  €  <SX^„  such  that 

(3.22)  ■®e(«(<j).fc)  (*/!>’’*»)  ~  V  Vh  € 

and 

(3.23)  =  0 

The  function  2)^  will  be  called  the  periodic  finite  element  solution.  Let  us  remark 
that  by  scaling  and  translating  the  mesh  we  can  solve  the  finite  element  problem 
in  the  master  cell  using  the  master  mesh. 

Lemma  5.  Let  (fi  €  h)^  and  let  ^  be  the  periodic  extension  of  tp  on 

S  Hi  <  El,  /  =  2,3.  Further  let  ijf  be  the  periodic  extension  of  on 

5(®®,ir/).  Then 


(3.24)  ff^^{ip,Vh)  V  Vfc  €  5J,o  ^/)) 

where 

(3.26)  5;^(S(.»,ft))  :=  |vfc  G  5fc(n)  :  Vh  has  its  support  in  5(»®,J?/)| 

Proof.  We  will  write 

(3.26)  A)  (**’’''*) 

(»j)e7i 


Using  the  periodicity  of  z^  and  ip  we  have 


(3.27) 


Z)  ^c(-(*.i).A)  (*A  > «a)  =  (^h .  Z  = 


=  Sci,o,H)  (^,  Z  (®  + 

(».i)€7^ 
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Realizing  that 


(3.28)  Xi  €  5J^.»(c(®®,h)) 

we  get  (3.24).  □ 

Let  us  now  use  9  =  p ,  where  p  has  been  defined  in  (3.15).  Then  by  lenuna  4 
we  have  tp  €  h)^  ^  p  =  p  end  denoting 

(3.29)  V’  :=  P  - 
we  have: 

Lemma  6. 

(3  JO)  (^,«»)  =  0  Vv^es^  (S(m‘,  Be)) 

and 

(3.3U)  C,  <  C,  ,  j  =  0,  1 

{3J16)  C,  A'*'-'  <  I IVV  1  <  C,  A-^'-^l/nAr ,  ;  =  0,  1 .  r  >  0 

where  Cg,  depend  on  px,  /13,  Hx  and  the  master  mesh. 

Proof.  We  have  (p  -  K^Vh)  =  ■Bsc-o.h,)  =  0  by 

lemma  5,  and  (3.30)  is  proven. 

From  (3.22)  and  lemma  4  we  get 

(3.32)  lll^lll,(.(y)ji)  =  III?  -  jJIIUjuia)  <  lll/>lll,(.l..)«  <  Ob'*' 
and  applying  the  duality  principle  we  get 

(3.33)  ll‘0lle(s(*>3),h)  —  ll*fc  ~  Plle(»(<.»),fc)  ^  ^^ll|p||le(«(**»),fc) 


15 


where  C  is  independent  of  p  and  h.  Realizing  that  there  are  Hj/h^  cells  in  5(0°,  Hi) 
we  get 


IllV'lllJc.o.ir,)  <  <  CHlh^ 


which  proves  the  right-hand  side  of  (3.31a). 

Because  =  C'{pi)  >  0  on  every  t  €  T^'^\  it  can  be  readily  seen  that 


from  which  the  left  hand  side  of  (3.31a)  follows. 

The  proof  for  the  norms  |  *  | ,  ||  *  j|  follows  analogous  steps  using  the  standard 
theory  of  X“  estimates.  Q 

Lemma  7.  Let  Q  be  a  polynomial  of  degree  (p  -f  1)  on  5(0°,  fTj),  fTj  <  Hi  and 
Qh^S\  ^5(0°,^])^  such  that 


(3.34a) 


inf 

X€5‘(5(.«',Jr,)) 


IIIQ  -  xlll5(»«,H,)  =  IIIQ  -  Qfc||ls(B*,iy,) 


(3.345)  /  (Q-  Qh)  =  0 

/5(a«  A) 

Further  let  H^  =  vH^,  0  <  tt  <  1,  dj  :=  if,  —  ^3  =  (1  —  'fr)H^.  Then 
(3-35)  WQ-Qh  =  II ) 

where  |A|  <  1  and  C  is  independent  of  h,  H-^  and  ilf. 

Proof.  We  have  Q  -  Qh-i^  =  +  —  Qh^  '5ji(5(0°,  ^2))  and 

Bs{^,h,){Q  -Qh-  V',Vfc)  =0  V  €  5J, 0(5(0°, ^2)) 
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and  ufling  Lemma  9  (see  below)  we  get 

II Q  -  Qfc  -  V'  B5(«*a)  ^  ^IIQ  -Qk-  V'll5(»»jr,)  ^ 

-  ^  0^^  “  ^kll5(«*A)  +  lt^lls(«»A))  -  ^^Ill^lll5(«»jsr,) 

The  above  follows  by  using  the  duality  principle  (5(0**,  if,)  is  convex  domain)  and 
by  observing  that 


lll<?  -  ^fclll5(«*4r,)  ^  lll^  “  -  IIIHIIs(»*4J,) 


and  that 

Moreover  |||p||ls(,o^j)  ^  ChPH^  and  thus 

(3.36) 

..  X.  ...  ^  ^ 

\\Q  Qh  V*  ll5(»».H,)  ^  ^  ji 

From  lemma  6 

II  ^  Il5(«®  A)  - 

and  hence 

HQ  Qh  V' lls{-o.H,) < c  11 V' lls(^.H,) 

Moreover 

II Q  ~  Qh  ll5(««.H,)  ^  II V*  ll5(*»,H,)  +  II Q  “  Qh  ~  ^  Use®", 

II Q  -  Qh  lls(®«,H,)  ^  II V'  lls(»«,H,)  ~  II Q  “  Qh  ~  lls(«» 

and  hence 


\\Q  -  Qh  lls(«'',Hj)  —  II  ^  lls(*o,Ji,)  ) 

which  completes  the  proof  of  the  lemma.  Q 
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We  will  now  summarize  the  theorems  on  interior  error  estimates  which  we  use 
in  this  Section.  For  the  proofs  see  [38]. 

Lemma  8.  Let  as  G  5(x®,iri),  :=  —  Bi  >  Ch.  Let 

(3.37«)  -«».»)=«  V  »  e  Hi.,)) 

Then 

ir.»(-)i  <  ^  I  “  -  I 

(3.376) 

Lemma  9.  Let  Wf^  G  5j(5(®®,fr<_i))  and 

(3.38a)  ~  ®  ^ €  5Xo(5(®®»^<_i)) 

Then 

(3.386)  ^  !l™'fcll5(«»,Hi_i) 

Let  us  now  prove  the  main  theorem  of  the  paper. 

Theorem  1 .  Let  <  H  <  H°  and  the  assumptions  I-lII  hold  with 

(3.39)  C,E?<K<C,H},  <.  =  5^.-  6(6;^.  />  =  ?+!-.  — 
Then  for  any  x  G  S{x^,H^) 

(3.40)  |^(,)|  =  |^(,)|  +  A«'-.. 

with  |A|  <  1  and  C  independent  of  h. 

Proof:  Let  JTj  =  (1  —  We  have 
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e*  =  «  -  ^  ^  +  Qfc  - 

and  hence 

II  e*  Bs(*»^,)  ^  II «fc  -  -  («  -  Q)  fi5(.*4i,)  +  H  Q  -  Qk  ls(«»4?,) 

Let  ns  estimate  ||  u,^  —  Q,^  —  (u  —  Q)  ||5(a«4r,)-  Noting  that 

-B5(.*  A)K  -  ^k  -  («  -  Q).  »)  =  0  V  V  €  52,0  (5(a»®,  U,)) 

and  using  lemma  8  we  get 

II  “k  -  ^k  -  -  Q)  Il5(.«,ir,)  <  j^(ll «  -  ^  -  X  ||s(.»^,)+ 

+  i  I  tt  -  Q  -  X 1 5(.«jr,))  +  ^IK  -  ^k  -  («  -  ^)ll5(^^,) 

Using  for  x  interpolant  of  u  —  Q  we  get 

II »  -  <3  -  X  lls(^ jr. )  <  Cff,*' ,  I «  -  «  -  X 1 5(^^, ,  < 

Further  we  have 

ll«k  -  Qk  -  ~  <?)lls(.»^,)  <  llttk  -  ■«»ll5(.»4r,)  +  llQk  -  Qll5(.0J5r,) 

and  using  assumption  III,  the  duality  principle  and  lemma  4  we  get  using  lemma  2 

IK  -<?»-(«-  «)ll5(„ji.)  < 

and  hence 

II «k -  ^k  ~  lls(««.ir,)  ^  C^HjhF  +  — — - h  — ^  ^ — J  <  Ch 

where 

H  =  mini  p+  —  ,  p+1  —  2*/,  ^  —  —  —  4i/,  p+1  —  —  —  4i/| 

la  a  a  ) 
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From  lemma  7  we  get 


Qh  il5(a*,Ha)  —  II  V*  |s(b»A)(^  + 


where  (Af  <  1.  Moreover  hram  lemma  6 


Hence 


II  «fc  Il5(-«4r,)  ^ 


and  analogously  the  other  inequality  so  we  get 


where 


II  Il5(a«^,)  —  II  V'  |l5(««,jsr,)(l  +  ACh^) 


u  =  mini  1  —  2i/,/9  —  —  —  4i/  —  p,  1  —  i  —  4i/| 
la  a  a  ) 


Assuming  that  /9  >  p  +  1  — 


0<e<<ro,  a=l  +  —  = 


1  6p+ 1 


6p  6p  ’ 


6(6p  + 1) 


we  get  with  Og  =  i/ , 


1  6p 

a  6p+ 1 


>  Og  ,  1  -  2|/  >  ffg  , 


B-1.-4V -p  =  v+l-e - - 7-^ — r  -  p  = 

^  a  ^  ^  6p+l  6(6p+l)  ^ 


36p  +  4  2  ^  1 

6(6p  4-1)  ^  6(6p  +  1)  ^  ~  6(6p  +1)  “  ’ 


=  _ 1 _ 

o  6p  +  1  6(6p  +  1) 


6(6p  +  1)  -  36p  -  4  2 

6(6p  +  1)  6(6p  +  1) 


when 
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We  thus  have 


1 

6(6p  +  l) 


=  e  =  i/  =  (r0>O. 


11  ®A  Ii5(««47«)  —  II  ^  +  ACh'*)  . 

Farther  we  have  for  any  m  6 

^(■)  =  ^((«  -  C)  -  K  -  %))(«)  +  -  «»)(*) 

and  so  proceeding  in  the  same  way  we  get 

□ 

Remark  3.6.  Let  us  now  comment  briefly  on  the  meaning  of  the  theorem.  We  will 
assume  that  the  mesh  is  periodic  in  a  small  re^on  in  the  interior  of  the  domain  and 
that  the  solution  is  smooth.  Outside  the  square  we  assume  neither  that  the  mesh 
is  periodic  nor  that  the  solution  is  smooth.  We  only  assume  that  the  pollution  in 
the  periodic  mesh-subdomain  is  controlled.  Then  the  superconvergence  is  directly 
related  to  the  function  ^  which  can  be  computed  in  the  master-cell. 

Remark  3.7.  Note  that  the  fact  that  we  have  used  square  cells  is  not  essential.  It 
is  only  essential  that  the  mesh  is  translation-invariant.  Let  us  also  underline  that 
the  assumption  that  K  is  the  identity  matrix  was  made  only  for  simplicity.  The 
theorem  holds  for  general  matrix  Ky  which  is  allowed  to  vary  smoothly  throughout 
the  domain. 

Remark  3.8.  Theorem  1  allows  us  to  make  a  computer-based  proof  of  the  existence 
of  the  superconvergence  points  for  the  case  when  F(u)  =  — —  (see  introduction). 

We  see  from  (3.40)  that  if  and  only  if  - — ^(»)  =  0  we  have 

axi 


(3.41) 


<^o  >  0 


and  C  is  independent  of  h.  Further  we  know  from  the  construction  of  r/f  that  (after 
scaling)  the  function  ‘0  is  independent  of  h.  Hence  the  relative  position  of  s;  in 


Q 

T  where  — — 0(a!)  =  0  is  independent  of  the  mesh. 
dXz 

exception  of  very  special  cases  we  have 


On  the  other  hand  with  the 
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(3.«) 

In  tliis  case  (which  is  always  assumed  in  the  classical  superconvergence  theories) 
onr  theorem  says  that  a  pomt  m  in  the  element  r  is  a  snperamveryence  point 

for  F{v)  —  if  and  only  if  -^(c)  =  0.  Because  we  are  interested  in  a  2/- 

OXj  OXf 

snperconvergence  point  for  U  being  a  sufficiently  large  set  of  solutions,  the  condition 
(3.42)  is  not  restrictive. 

Remark  S.9.  We  discussed  above  only  the  heat-conduction  problem.  We  will 
address  the  superconvergence  for  the  dastidty  problem  in  a  forthcoming  paper 
[30]. 

RemaH:  S.IO.  Although  in  general  the  superconvergence  points  do  not  exist,  a  re¬ 
gion  of  points  where  the  accuracy  of  the  gradient  is  much  larger  than  the  maximal 
error  in  the  element  still  could  exist.  This  question  is  addressed  in  forthcoming  pa¬ 
pers  ([28],  [29])  in  which  we  introduce  and  analjrxe  the  notion  of  q-snperconvergence 
points,  0  <  q  <  1.  The  case  q  =  0  coincides  with  the  classical  superconvergence 
(see  also  (1.4)). 

4  The  methodology  for  finding  the  supercon¬ 
vergence  points 

In  this  paper  superconvergence  is  treated  as  local  behavior  and  hence  we 
are  making  assumptions  about  the  local  behavior  of  the  solution  in  the  interior  of 
the  domain.  We  will  consider  the  class  of  solutions  which  are  locally  smooth  in 
5(si*’,  H)^  namely, 

(4.1)  :=  {«  €  jr'(n)|  I  D-«  1 <  Jr ,  0<|a|<p+2| 

where  5(k°,J7)  denotes  an  interior  subdomain  of  interest  in  which  the  mesh  is 
locally  periodic  as  described  above  (the  subdomain  must  be  a  finite  distance  away 
from  the  boundary  and  points  of  roughness  of  the  source  term).  In  many  instances 
we  are  only  interested  in  the  subclass  of  solutions  in  which  are  "harmonic” 
(we  say  that  u  is  "harmonic”  if  it  satisfies  the  homogeneous  differential  equation), 
namely. 


(4.2) 


:=  e  W®  I  C{u)  =  0  in  flj 
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We  may  also  assume  that  the  functions  are  "harmonic”  in  a  subdomain  which  is 
slightly  bigger  than  S(x^,H)  and  which  includes  S{m^,H)  in  its  interior. 

In  the  previous  Section  we  proved  that  we  can  obtain  the  asymptotic  values 
of  the  error  for  any  smooth  solution  in  the  interior  of  a  periodic  mesh-subdomain 
by  solving  a  periodic  boundary-value  problem,  using  the  master-mesh  T  over  the 
master-cell  c,  with  data  obtained  from  the  local  (p  -i-  l)-degree  Taylor-series  ex¬ 
pansion  of  the  exact  solution.  Based  on  this  result  we  will  construct  a  numerical 
procedure  to  determine  the  superconvergence  points  (which  may  not  exist)  for  a 
given  class  of  smooth  solutions  by  employing  the  corresponding  class  of  (p  -|-  1)- 
degree  monomials. 

a.  The  classes  of(p  +  l)-degree  monomial  solutions. 

Let  us  assume  that  for  a  given  locally  periodic  grid  with  corresponding  periodic 
master-mesh  T,  given  material  orthotropy  and  given  class  of  smooth  solutions  li 
we  consider 

(4.3)  Q  :=  |t?|  Q{x^^x^  =  51  ~  } 

'  k=i  t=o  * 

the  class  of  (p-h  l)-degree  monomials  which  occur  in  all  (p -I-  l)-degree  Taylor-series 
expansions  about  a:*’  of  functions  from  U.  Here  h  =  1, . . .  ,nd  denotes  a  set  of 
linearly  independent  monomials  which  span  Q.  For  example  let  us  assume  that  U 
is  the  class  of  smooth  solutions  U°  given  in  (4.1);  in  this  case  we  may  choose 

(4.4)  *a)  •=  ,  l<fc<nd  =  p-|-2 

and  we  obtain  the  class  of  all  (p  -I-  l)-degree  monomials  Q?. 

In  many  practical  applications  we  are  interested  only  in  the  class  of  "harmonic” 
solutions  11“^”  given  in  (4.2).  Then  Q  is  the  two-dimensional  space  of  “harmonic” 
monomials  and  will  be  denoted  by  Q"^.  For  example  in  the  case  of  Laplace’s 
equation  =  =  1,2)  we  have 

(4.5a)  :=  |q^|  =  53  Q? (®i>®2)  }  > 

*=i  ^ 

(4.55)  Qf  (xi,®,)  =  ne{z^^) ,  Qf  (xi,®j)  =  ,  z  =  z,+ix,. 

Note  that  in  this  case  we  do  not  use  the  quotes  since  the  functions  are  harmonic 
in  the  classical  sense  (i.e.  they  satisfy  Laplace’s  equation). 
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b.  The  error  in  the  periodic  finite-element  solution. 

For  a  given  exact  solution  Q  €  Qy  material-orthotropy  matrix  Ky  periodic 
grid  T  we  obtain  the  error  in  T  by  solving  the  following  periodic  boimdary-value 
problem: 

Find  €  H^{c)  such  that 

(4.6a)  BtizJiyV^)  =  Bi(pyv^)  V  6  ^.pbr(c) 

where  p,  •SXj>br(c)  are  defined  in  (3.15),  (3.20),  respectively,  and 

(4.66)  = 

and  we  obtun  the  error 


(4.7)  p-z^. 

•  Let  F(u)  denote  the  solution-quantity  of  interest;  for  example  below  we  will 

.1  \  du  V  /X  du  du 

consider  F(n)  =  ^  or  F(u)  =  ^^(it)  = 

c.  Superconvergence  points  of  F{u)  for  the  class  Q. 

We  will  say  that  s  is  a  superconvergence  point  in  the  element  r  6  T  for  the 
functional  F{u)  and  the  class  of  solutions  Q  if 

(4.8)  =  0  V  Q  e  C 

Here  rl>  is  the  error  which  corresponds  to  the  solution  Q. 

To  determine  the  superconvergence  points  for  F(tt),  in  r  €  T  for  the  class  Q  we 
define  the  zero-contour  for  F(V’*)  corresponding  to  the  linearly  independent  mono¬ 
mials  Qi ,  i  =  1 . . .  ,nd  in  the  element  r  6  T  and  for  the  material-orthotropy  K 

(4.9)  ;  T,  f ;  JT)  ;=  {  *  €  r  I  F{r)  (»)  =  0  }  , 

where  is  a  basis  for  Q  and  i/)'  denotes  the  error  corresponding  to  Q^.  The 

set  of  superconvergence  points  in  the  element  r  is  given  by 
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(4.10) 


i=l 

In  the  implementations  the  contours  r,  T;  K)  were  approximated  by 

piecewise  linear  continuous  lines  in  each  element  r  using  a  background  mesh  of 
linear  triangles.  The  initial  approximations  of  the  superconvergence  points  were 
obtained  from  the  intersections  of  the  approximate  contours.  In  order  to  find  more 
accurate  values  of  the  superconvergence  points  we  constructed  the  function 

(4.11)  G„.,(42;T.f;  JST)  :=  ^ 

t=l 


The  superconvergence  points,  if  they  exist,  are  located  at  the  zeros  of  this  function. 
A  minimization  algorithm  was  employed  to  improve  the  initial  approximations  of 
the  superconvergence  points  obtained  from  the  intersections  of  the  approximate 
contours  by  minimizing  ^F(uy 

Remark  Note  that  the  superconvergence  points  for  the  functional  F{y)  and 
the  class  Q  may  not  exist  in  the  element  r  i.e.  the  set  Xp^^{Q\TyT\K)  may 
be  empty.  This  is  indeed  the  case  when  the  function  aoes  not  have  any 

zeros.  In  the  case  that  Q  is  the  class  of  harmonic  monomials,  we  have  only  two 
families  of  zero-contours  which,  in  general,  will  intersect.  Thus  we  can  always  find 
superconvergence  points  for  the  class  of  harmonic  solutions. 


5  Numerical  study  of  superconvergence  for  pe¬ 
riodic  meshes  of  triangles  and  squares 

We  will  now  demonstrate  that,  by  using  the  methodology  outlined  in  the 
previous  Section,  we  can  find  the  superconvergence  points  for  the  components  of 
the  gradient  (or  the  components  of  the  flux  in  the  case  of  orthotropic  material)  in 
the  interior  of  periodic  meshes  of  triangular  or  square  elements.  In  the  niunerical 
examples  we  will  study  the  following  questions: 

1.  For  periodic  meshes  of  triangles  with  various  mesh- topologies  and  elements 
of  degree  p  where  are  the  superconvergence  points  for 

a.  the  class  of  harmonic  smooth  solutions; 

b.  the  class  of  general  smooth  solutions? 
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2.  For  meshes  of  square  elements  of  degree  p  where  are  the  superconvergence 
points  for 

a.  the  tensor-product  space 

b.  the  serendipity  space  ^(r); 

c.  the  intermediate  space  S*^(f)? 

We  will  answer  these  questions  using  the  numerical  approach. 

Remark  5.1.  We  remark  that  for  unsmooth  solutions  the  superconvergence  points 
do  not  exist,  in  general. 

Remark  5.2.  Although  in  Section  4  we  addressed  only  the  case  of  triangular  ele¬ 
ments  the  theory  holds  for  the  general  element-spaces  we  address  below. 

5.1  Determination  of  the  superconvergence  points  for  pe¬ 
riodic  meshes  of  triangles 

The  majority  of  the  results  for  the  superconvergence  points  for  the  triangular 
elements  in  the  literature  are  given  exclusively  for  the  Regular  pattern  (which  is 
also  known  as  the  ihrtt~dirtciional  mesh  and  is  shown  in  Fig.  2a)  and  linear  and 
quadratic  elements;  see  for  example  [11],  [13],  [15-26].  To  our  knowledge  there  are 
no  systematic  results  reported  in  the  literature  about  the  superconvergence  points 
for  other  mesh-pattems  and  p  >  3.  Here  we  will  employ  the  numerical  approach  of 
Section  4  to  find  the  superconvergence  points  for  all  the  patterns  shown  in  Fig.  2 
and  for  elements  of  degree  p,  1  <  p  <  7.  Moreover  we  will  find  the  superconvergence 
points  for  the  harmonic  class  and  the  general  class  of  smooth  solutions. 

We  determined  the  superconvergence  points  for  the  -derivative  of  the  gradient 
for  the  class  of  harmonic  solutions  by  locating  the  points  of  intersection  of  the 
contours  C*.  ;r,T;  J),  t  =  1,2.  We  then  checked  if  these  points  are  also 

superconvergence  points  for  the  class  of  general  solutions  by  computing  the  values 
of  the  Xj-derivative  of  the  error,  — ,  corresponding  to  the  linearly  independent 

UX\ 

monomials  which  are  needed  to  complete  the  basis  of  the  harmonic  space  Q,^  into 
the  basis  of  the  general  space  Q^. 

In  Table  1  we  give  the  set  of  superconvergence  points  X  ^  T;  I)  for  the 

- 

class  of  harmonic  solutions  and  for  the  periodic  meshes  T  shown  in  Figs.  2a,  2b,  2c 
and  the  element  r  shown  with  gray-shading  in  these  Figures  (note  that  the  points 
are  the  same  for  all  three  patterns).  The  superconvergence  points  for  the  class  of 
harmonic  solutions  for  the  element  (resp.  r,  )  in  the  Criss-Cross  pattern  shown 
in  Fig.  2d  are  given  in  Table  2  (resp.  Table  3).  In  Tables  4a-4d  (resp.  5a-5b)  we 
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indicate  which  points  for  Table  1  (resp.  Tables  2,  3)  are  also  superconvergence 
points  for  the  class  of  general  solutions  . 

In  Figs.  3a-3d  we  give  a  schematic  example  of  how  the  superconvergence  points 
for  the  class  of  harmonic  solutions  in  the  cubic  elements  and  the  four  mesh-pattems 

are  obtained  from  the  intersections  of  the  zero-contours  for  —  for  the  two  linearly- 

OX\ 

independent  harmonic  monomials  of  degree  (p  -|- 1).  Figs.  3a-3d  also  show  how  to 
interpret  the  points  given  in  Tables  1-3  for  all  the  elements  in  each  pattern.  In 
Fig.  4a  we  show  the  contours  of  zero-error  in  the  x^-derivative  for  the  two  harmonic 
monomials  for  the  shaded  element  from  the  Regular-pattern  while  in  Fig.  4b  we 
also  included  the  zero-contours  for  the  monomials  which  are  needed  to  complete 
the  harmonic  basis  into  the  general  one.  We  note  that  there  is  only  one  point  which 
is  common  to  the  zero-contours  of  all  the  monomials,  namely  the  midpoint  of  the 
side  of  the  triangle  which  is  parallel  to  the  x^-axis.  As  indicated  in  Table  4a,  this 
point  is  the  only  superconvergence  point  (for  the  class  of  general  solutions)  for  the 
Xj -derivative  in  element  r  in  the  Regular  pattern. 

We  make  the  following  observations: 

a.  For  the  harmonic  solutions  there  are  always  superconvergence  points  for  the 
Xj-derivative  for  all  the  mesh-pattems.  The  number  of  the  points  increases 
with  the  degree  p  of  the  elements. 

b.  Only  very  few  of  the  superconvergence  points  for  the  class  of  harmonic  solu¬ 
tions  are  also  superconvergence  points  for  the  class  of  general  solutions.  For 
example 

(i)  For  even  p’s,  p  >  4  there  are  no  superconvergence  points  for  the  class  of 
general  solutions. 

(ii)  For  the  odd  p’s,  p  >  3  and  the  Regular  and  Chevron  pattern  there  is 
only  one  superconvergence  point  for  the  class  of  general  solutions.  It  is 
located  at  the  mid-point  of  the  side  parallel  to  the  Xj-axis. 

(iii)  For  the  Union-Jack  pattern  there  are  no  superconvergence  points  for  the 
class  of  general  solutions  for  elements  of  any  degree  p. 

c.  The  superconvergence  points  for  the  x,- component  of  the  gradient  for  the 
mesh-patterns  shown  in  Fig.  2  can  be  obtained  using  synunetry  arguments. 

We  now  demonstrate  that  when  the  points  given  in  Tables  1-3  are  employed  in 
model  computations  the  superconvergence  of  the  x^-derivative  can  be  observed  for 
relatively  coarse  meshes.  To  study  the  superconvergence  in  actual  computations 
we  employed  the  discrete-norm  suggested  in  [8,9],  namely 
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Here  denotes  the  superconvergence  points  in  the  triangles  of  the  mesh 

which  participate  in  the  definition  of  the  discrete  norm. 

We  considered  the  domun  fl  =  (0|1)^  which  is  meshed  by  a  uniform-grid  of 
elements  in  the  Regular  pattern.  We  computed  the  values  of  the  discrete  norm  us¬ 
ing  points  from  the  elements  in  the  subdomain  Qo  =  (0.25,0.75)’.  We  considered 

meshes  with  cell-size  h  =  we  computed  three  types  of  discrete 
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norms  E',  E'".  In  particular  E'  includes  the  superconvergence  points  given  in 

Table  1  for  each  element  in  fig,  E"  includes  only  those  points  for  each  element  in  fig 
which  are  superconvergence  points  for  the  Xj -derivative  for  the  class  of  general  so¬ 
lutions  (see  Table  4)  and  E‘"  includes  the  points  from  the  elements  in  fig  which  are 
superconvergence  points  only  for  the  class  of  harmonic  solutions.  Thus  we  expect 
to  see  that  all  three  norms  are  superconvergent  if  the  exact  solution  is  harmonic 
while  only  E"  is  also  superconvergent  for  solutions  which  are  not  harmonic. 

In  Table  6a  (resp.  Table  6c)  we  give  the  values  of  h~^E',  h~^E",  h~^E'" 
(resp.  h~*E*t  h~*E",  h~*E"')  computed  from  the  solution  of  a  Dirichlet  boundary- 
value  problem  in  fl  using  quadratic  (resp.  cubic  elements)  with  data  consistent 
with  the  harmonic  solution  u^{x^^x^)  =  sin(irxi)sinh(xx2).  We  note  that  as  h 
is  subdivided  the  values  h-^E\  h’^E",  h-^E‘"  (resp.  h-*E',  h'^E",  h-*E'‘')  con¬ 
verge  to  constants;  this  shows  the  superconvergence  of  all  three  versions  of  the 
discrete-norm. 

In  Table  6b  (resp.  Table  6d)  we  give  the  values  of  h~^E\  h~^E",  h~^E'" 
(resp.  h~^E\  h~*E‘\  h~^E"')  computed  from  the  solution  of  a  Dirichlet  boundary- 
value  problem  in  fl  using  quadratic  (resp.  cubic)  elements  with  data  consistent  with 
the  solution  u^(zj,x,)  =  sin(TX^)sin(wx2)  (which  is  not  harmonic).  We  note  that 
as  h  is  subdivided  the  values  of  h~^E\  k~^E",  h~^E"'  (resp.  h~^E\  h~*E",  h~^E'") 
converge  to  constants;  this  shows  that  only  the  values  of  E"  are  superconvergent 
in  this  case. 

5.2  Determination  of  the  superconvergence  points  for  pe¬ 
riodic  meshes  of  squares 

We  also  employed  the  computer  based  approach  to  find  the  superconver¬ 
gence  points  for  meshes  of  square  elements.  In  particular  we  studied  the  effect 
of  the  choice  of  the  finite-element  space  (tensor-product,  serendipity,  intermedi¬ 
ate  spaces  as  defined  in  (2.9),  (2.10),  (2.11))  on  the  superconvergence  points  for 
the  x^-derivative.  For  the  tensor-product  space  it  is  well  known  that  the  points 
of  the  (p  X  p)  Gauss-Legendre  product  rule  are  the  superconvergence  points  for 
the  x^-  and  X2-derivative  simultaneously  for  elements  of  degree  p  (see  for  example 
[9]).  To  our  knowledge,  superconvergence  points  for  the  serendipity  family  and  the 
intermediate  have  not  been  reported  in  the  literature  for  p  >  3. 
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We  used  the  numerical  approach  to  determine  the  superconvergence  points  for 
the  sB^-derivative  for  the  case  of  Poisson’s  equation.  We  observed  the  following: 

a.  For  the  tensor-product  and  the  intermediate  family  of  element-spaces,  su¬ 
perconvergence  of  the  Zj-derivative  occurs  along  the  constant  z^-Gauss  lines 
which  pass  through  the  Gauss- Legendre  points  of  degree  p  along  the  x^-axis. 

b.  For  the  serendipity  family  the  superconvergence  points  for  the  x^-derivative 
are  given  as  follows: 

(i)  For  p  =  1,2  and  the  class  of  general  solutions  the  points  are  exactly  the 
same  Gauss-lines  as  for  the  tensor-product  and  the  intermediate  family. 

(ii)  For  p  =  3  and  the  class  of  general  solutions  four  superconvergence  points 
and  one  superconvergence  line  (the  line  Xy  =  0)  were  found,  as  reported 
in  Table  7  and  shown  in  Fig.  5.  For  the  class  of  harmonic  solutions  three 
superconvergence  lines,  which  are  shown  with  thick  continuous  line  in 
Fig.  5,  were  obtained. 

(iii)  For  p  >  4  there  are  no  superconvergence  points  for  the  class  of  general 
solutions.  For  the  class  of  harmonic  solutions  several  superconvergence 
points  exist;  these  points  are  listed  in  Table  7  and  shown  schematically 
in  Figs.  6a-6d  for  p  =4,  5,  6,  7,  respectively. 

To  demonstrate  the  superconvergence  of  the  Xj-derivative  for  the  points  in 
Table  7  for  p  =  3  and  4,  we  considered  a  Neumann  boundary- value  problem  in 
SI  =  (0, 1)’  with  data  consistent  with  the  harmonic  solution  and  the  general 
solution  u^.  We  computed  the  discrete  norm  E  using  points  &om  the  elements  In 
the  subdomain  fig  =  (0.25,0.75)^. 

In  Table  8a  we  give  the  value  of  h~*E  where  the  (3  x  3)  Gauss- Legendre  points 
from  each  element  in  Qq  were  used  in  the  definition  of  the  discrete  norm  E]  the 
Neumann  problem  with  data  obtained  from  the  general  solution  was  solved 
with  elements  of  degree  p  =  3  using  the  tensor-product  space  and  the  intermediate 
space.  We  observe  that  for  both  spaces  the  values  of  E  are  superconvergent,  as 
expected. 

In  Table  8b  we  give  the  values  of  E  and  h~^E  (resp.  h~*E)  computed  from 
the  finite  el^ement  solution  of  the  Neumann  problem  with  data  consistent  with  the 
general  sol.  tion  and  using  the  (3  x  3)  Gauss-Legendre  points  from  each  element 
in  Go  (resp.  the  four  points  given  in  Table  7  and  the  point  (0,0)).  We  note  that  the 
values  of  E  are  superconvergent  only  when  the  points  from  Table  7  are  employed. 

In  Table  8c  we  give  the  values  of  h~’^E  (resp.  h~*E)  computed  from  the  finite 
element  solution  of  the  Neumann  problem  with  data  consistent  with  the  harmonic 
solution  (resp.  the  general  solution  u^).  Here  we  used  quartic  serendipity  finite 
element  space  and  the  points  given  in  Table  7  for  p  =  4.  We  note  that  the  values 
of  the  discrete  norm  E  are  superconvergent  only  in  the  case  of  the  harmonic. 
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5.3  Effect  of  the  material  orthotropy  on  the  superconver¬ 
gence  points  for  the  flux 

In  the  above  examples  we  gave  the  snperconvergence  points  for  the  components 
of  the  derivative  of  finite  element  solution  for  Poisson’s  equation.  We  now  study 
the  effect  of  orthotropy  on  the  superconvexgence  points. 

We  considered  square  cubic  elements  (p  =  3).  We  let  the  material-orthotropy 

matrix  K  =  K{6)  where  K{6)  has  ratio  of  principal  values  =  10  and  6  de- 

notes  the  angle  (measured  counterclockwise)  between  the  z^-axis  and  the  principal- 
axis  of  orthotropy  corresponding  to  Kaia  .  We  considered  the  class  of  “harmonic” 
solutions  and  the  corresponding  space  Q“^“  of  functions  which  can  be  ex¬ 
pressed  as  linear  combinations  of  the  “harmonic”  monomials 


(5.2)  x”  =  x;  +  <x;. 


Here  denotes  the  axes  of  the  coordinate-system  for  which  the  differential 

operator  is  transformed  to  isotropic  i.e. 


(5.4a) 

where 


(5.46)  v!\xl,xl)  :=  ^(xijXj) . 

In  Fig.  7  we  give  the  zero-contours  C',^(„)(^<^”;T,T';fir),  i  =  1,  2,  and  the 

superconvergence  points  for  the  tensor-product  space  and  the  intermediate 
space  S'^.  We  note  that  the  zero-contours  for  the  two-spaces  are  identical  and  the 
superconvergence  points  do  not  depend  on  the  material-orthotropy.  (The  super¬ 
convergence  points  for  the  components  of  g(u)  for  elements  of  degree  p  from  the 
tensor-product  or  the  intermediate  space  are  located  at  the  (p  x  p)  Gauss- Legendre 
points  in  the  element.) 

In  Fig.  8  we  give  the  zero-contours  ■*^)>  *  =  1»  2,  and  the  super- 

convergence  points  for  the  cubic  serendipity  element  and  various  grid-material  ori¬ 
entations.  We  note  that,  in  general,  the  points  depend  on  the  material-orthotropy 

(on  the  ratio  and  the  grid-material  orientation  6).  This  dependence  should 
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also  be  expected  for  serendipity  elements  of  higher-degree  (p  >  4). 
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In  Figs.  9a-9h  we  give  the  zero-contours  r,  T;  K),  i  =  1,  2,  and  su¬ 

perconvergence  points  for  quadratic  triangular  elements  and  various  grid-material 
orientations  [  =  10,  €  [0“,90“]V  In  Figs.  9a-9h  we  show  the  zero  contours 

and  superconvergence  points  for  the  element  r  (shaded  gray  in  Fig.  2a)  for  the 
Regular  pattern.  From  the  Figures  it  is  clear  that  the  location  and  number  of 
superconvergence  points  depends  on  the  material-orthotropy.  Note  that  for  the 
given  orthotropy  and  0  =  0°  there  are  two  points  of  superconvergence;  for  6  =  90° 
there  are  four  superconvergence  points  while  for  any  of  the  grid-material  orienta¬ 
tion  there  is  only  one  superconvergence  point.  Based  on  the  above  results  we  can 
conjecture  that  for  meshes  of  triangles  the  number  and  location  of  the  supercon¬ 
vergence  points  for  the  ‘‘harmonic”  class  of  solutions  depend  on  the  grid-material 
orientation. 


6  Summary  of  conclusions. 

1.  We  presented  a  computer-based  methodology  which  can  be  employed  to  de¬ 
termine  the  superconvergence  points  for  the  derivatives  in  the  interior  of  any 
locally  periodic  mesh  of  triangles  or  squares  of  any  degree  p. 

2.  The  methodology  takes  directly  into  account  the  topology  of  the  mesh,  the 
element  polynomial  spaces  and  the  nature  of  the  solution. 

3.  We  employed  the  methodology  to  determine  the  superconvergence  points  for 
several  cases  of  practical  interest.  We  made  the  following  observations: 

a.  For  the  class  of  harmonic  solutions  the  superconvergence  points  always 
exist  for  any  mesh-pattem  and  type  of  elements. 

b.  For  the  class  of  general  solutions  the  superconvergence  points  may  not 
exist. 

c.  For  the  square  elements  and  tensor-product  and  the  intermediate  element- 
spaces  of  degree  p  the  superconvergence  points  are  the  points  of  the 
(p  X  p)  Gauss-Legendre  product  rule. 

d.  For  square  elements  and  the  serendipity  space  of  degree  p  the  super¬ 
convergence  points  depend  on  the  class  of  solutions  and  the  material- 
orthotropy  for  elements  of  degree  p  >  3. 

e.  For  triangular  elements  the  superconvergence  points  depend  on  the 
mesh-pattem,  the  class  of  solutions  and  the  material-orthotropy. 
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List  of  Figures 

Figure  1.  Locally  periodic  grids,  (a)  An  example  of  a  grid  with  a  locally  periodic 
subdomain  meshed  with  the  Regulat  pattern  shown  in  Fig.  2a;  (b)  The  subdomain 
5(89°,  J7)  with  the  subdomains  5(0**,  J7^),  /  =  1,2,  in  its  interior. 

Figure  2.  Periodic  meshes  of  triangles:  (a)  Regular  pattern;  (b)  Chevron  pattern; 
(c)  Union-Jack  pattern;  (d)  Criss-Cross  pattern. 

Figure  3.  Superconvergence  points  for  the  x^-derivative  for  the  class  of  harmonic 
solutions  of  Poisson’s  equation:  Cubic  triangular  elements.  The  superconvergence 
points  are  located  at  the  intersection  of  the  zero-contours  C  am  (Qf;  r,  T;  /),  t  = 

1,2,  for  the  x^-derivative  of  the  error  for  the  harmonic  monomial  solutions.  In 
the  Figures  C ^  (Qf:r^f:I)  (resp.  C ^  (Qf :  r.  f\  I))  was  drawn  using  thin  (resp. 

thick)  continuous  line  in  each  element  r  (the  element  boundaries  were  drawn  using 
dashed  line).  The  zero-contours  and  the  superconvergence  points  are  given  for:  (a) 
The  Regular  pattern;  (b)  The  Chevron  pattern;  (c)  The  Union-Jack  pattern;  (d) 
The  Criss-Cross  pattern. 

Figure  4.  Superconvergence  points  for  the  x^-derivative  for  the  classes  of  harmonic 
and  general  solutions  of  Poisson’s  equation:  Cubic  triangular  element  r  from  the 
Regular  pattern,  (a)  The  triangle  with  the  zero-contours  C^(Qf  ;t,T;  J),  i  = 

1,2  for  the  quartic  harmonic  monomiab;  (b)  The  triangle  with  the  zero  contours 
TyT] I),  i  —  1, . . . , 4  for  all  the  quartic  monomials  for  which  the  error  in 
the  x^-derivative  does  not  vanish  identically.  The  zero  contours  for  the  harmonic 
monomials  where  drawn  using  thick  continuous  line;  the  zero-contours  for  the  other 
two  monomials  were  drawn  using  thin  continuous  line.  Note  that  the  midpoint 
of  the  horizontal  side  of  the  triangle  belongs  to  the  zero-contours  for  all  quartic 
monomials  and  is  a  superconvergence  point  (the  only  one)  for  general  solutions. 

Figure  5.  Cubic  serendipity  element:  Superconvergence  points  for  the  x^-derivative 
for  the  class  of  general  solutions  of  Poisson’s  equation.  The  superconvergence  points 
are  located  at  the  intersection  of  the  zero-contours  (Q^;  r,  T;  J),  t  =  1,2,  where 

Qj(xj,X3)  =  x^  and  Q,(X},x,)  =  x^  —  6xfx|  -|-  x^.  The  zero-contours  for 
(resp.  Q2)  were  drawn  with  thin  (resp.  thick)  continuous  line.  Note  that  there 
are  four  distinct  superconvergence  points  and  a  superconvergence  line  (the  line 
indicated  by  a  multitude  of  points  and  passes  through  the  center  of  the  element). 

Figure  6.  Serendipity  elements  of  degree  p  (4  <  p  <  7)  :  Superconvergence 
points  for  the  x^-derivativefor  the  class  of  harmonic  solutions  of  Poisson’s  equation. 
The  superconvergence  points  are  located  at  the  intersection  of  the  zero-contours 
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(Qf ;  Tyf]I)i  =  1,2,  for  the  Xj-deriTative  of  the  error  for  the  harmonic  mono¬ 
mial  solutions.  In  the  Figures  (resp.  C^{Q^\T,f\I)  was  drawn 

using  thin  (resp.  thick)  continuous  line.  Zero-contours  and  superconvergence 
points  for  the  x^ -derivative  for  the  serendipity  square  element  of  degree:  (a)  Four; 
(b)  Five;  (c)  Six;  (d)  Seven.  Note  that  there  are  no  superconvergence  points  for 
the  class  of  general  solutions  for  these  dements. 

Figure  7.  Effect  of  the  material  orthotropy.  Superconvergence  points  for  the 

Xj -component  of  the  flux,  qi(tt)  =  - h  /fijo —  for  the  class  of  “harmonic” 

solutions  of  the  equation  of  orthotropic  heat-conduction:  Cubic  elements.  The 
superconvergence  points  are  located  at  the  intersection  of  the  zero-contours 
C,  („)(QJ®";t,T;  JiC(d)),  t  =  1,2.  In  the  Figures  the  contours  were  drawn  for  ma- 
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trices  K(6)  with 


=  10  and  for  grid-material  orientations  6  =  0°,  15'*,30'’,45'’, 


60°,  75°,  90°.  The  grid-material  orientation  is  measured  by  the  angle  6  between 
the  x^-axis  and  the  principal-axis  corresponding  to  (measured  counterclock¬ 
wise).  (a)  Zero-contours  and  superconvngence  points  for  the  tensor-product  fam¬ 
ily;  (b)  Zero-contours  and  superconvergence  points  for  the  intermediate  family. 
Note  that  the  two  Figures  are  identical.  The  locations  of  the  superconvergence 
points  (which  are  at  the  (3  x  3)-Gauss-Legendre  points)  do  not  depend  on  the 
grid-material  orientation. 

Figure  8.  Effect  of  the  material  orthotropy.  Superconvergence  points  for  the  x^- 
component  of  the  flux,  gj(«)  =  - h  —  for  cl“s  of  “harmonic”  solu- 

vXj  vXj 

tions  of  the  equation  of  orthotropic  heat-conduction:  Cubic  serendipity  elements. 
The  superconvergence  points  are  located  at  the  intersection  of  the  zero-contours 
G,  (u)(Qi^”‘i‘r,r;  Jir(d)),  t  =  1,2.  In  the  Figures  the  contours  were  dravm  for 

matrices  K{6)  with  -*°”  =  10  and  for  grid-material  orientations  6  =  0.  The 

■/V  Twivt 

grid-material  orientation  is  measured  by  the  angle  6  between  the  Xj-axis  and  the 
principal-axis  corresponding  to  Kjaia  (measured  counterclockwise).  Zero-contours 
and  superconvergence  points  for  grid  material  orientation  6:  (a)  6  =  0°;  (b)  6  =  15°; 
(c)  0  =  30°;  (d)  6  =  45°;  (e)  6  =  60°;  (f)  6  =  75°;  (g)  0  =  90°.  Note  that  for 
6  =  0°,90°  there  are  three  superconvergence  lines.  The  zero- contours  and  the 
superconvergence  points  for  all  the  above  0h  are  displayed  in  Fig.  8h.  Note  that 
there  is  only  one  superconvergence  point,  located  at  the  center  of  the  element, 
which  is  valid  for  all  grid-material  orientations. 

Figure  9.  Effect  of  the  material  orthotropy.  Superconvergence  points  for  the 

du  du 

Xj-component  of  the  flux,  ^^(u)  =  K^— - h  K^2'a —  fo^  “harmonic” 
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solutions  of  the  equation  of  orthotropic  heat-conduction:  Quadratic  triangular 
elements,  Regular  pattern.  The  superconvergence  points  are  located  at  the  inter¬ 
section  of  all  the  zero-contours  i  =  1,2.  In  the  Figures 


the  contours  were  drawn  for  matrices  K{0)  with 


=  10  and  for  grid-material 


orientations  0,  where  9  is  the  angle  between  the  x^-axis  and  the  principal-axis  cor¬ 
responding  to  KjBoa  (measured  counterclockwise)  for  the  element  r  (shaded  gray  in 
Fig.  2a).  Zero-contours  and  superconvergence  points  for  (a)  6  =  0°;  (b)  9  —  15°; 
(c)  9  =  30°;  (d)  9  =  45°;  (e)  9  =  60°;  (f)  9  =  75°;  (g)  9  =  90°.  Note  that 
for  0  =  0°  there  are  two  superconvergence  points;  the  9  =  90°  there  are  four  su¬ 
perconvergence  points.  The  zero- contours  and  superconvergence  points  for  all  the 
angles  9,  mentioned  above,  are  plotted  in  Fig.  9h.  Note  that  there  is  no  common 
superconvergence  points  for  all  the  grid-material  orientations. 
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List  of  Tables 


Table  la.  Superconvergence  points  for  periodic  meshes  of  triangular  elements 
of  degree  p  (1  <  p  <  S).  Harmonic  solutions:  Regular,  Chevron  and  Union- 
Jack  patterns.  The  superconvergence  points  for  the  Xi-derivative  are  listed  for  the 
shaded-element  showm  in  Figs.  2a,  2b,  2c.  Note  that  the  superconvergence  points 
in  the  element  are  the  same  for  all  three  patterns. 

Table  lb.  Superconvergence  points  for  periodic  meshes  of  triangular  elements 
of  degree  p  (6  <  p  <  7).  Harmonic  solutions:  Regular,  Chevron  and  Union- 
Jack  patterns.  The  superconvergence  points  for  the  Zi-derivative  are  listed  for  the 
shaded-element  shown  in  Figs.  2a,  2b,  2c.  Note  that  the  superconvergence  points 
in  the  element  are  the  same  for  a'',  three  patterns. 

Table  2a.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  (1  ^  p  ^  5).  Harmonic  solutions:  Criss-Cross  pattern.  The  supercon¬ 
vergence  points  for  the  xi -derivative  are  listed  for  the  shaded-element  ri  shown  in 
Fig.  5d. 

Table  2b.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  =  6.  The  superconvergence  points  for  the  Xi-derivative  are  listed  for  the 
shaded-element  ti  shown  in  Fig.  2d. 

Table  3a.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  (1  <  p  <  5).  Harmonic  solutions:  Criss-Cross  pattern.  The  supercon¬ 
vergence  points  for  the  zi -derivative  are  listed  for  the  shaded-element  rj  shown  in 
Fig.  2d. 

Table  3b.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  =  6.  The  superconvergence  points  for  the  zi-derivative  are  listed  for  the 
shaded-element  rj  shown  in  Fig.  2d. 

Table  4a.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p(l^p<7).  General  solutions:  Regular  pattern.  This  Table  refers  to  the 
superconvergence  points  given  in  Tables  la,  lb  for  the  class  of  harmonic  solutions. 
The  points  of  Tables  la,  lb  which  are  also  superconvergence  points  for  general 
solutions  for  this  pattern  are  indicated  by  the  word  yes. 

Table  4b.  Superconvergence  points  for  periodic  meshes  of  triang\ilar  elements  of 
degree  p  (1  <  p  <  7).  General  solutions:  Chevron  pattern.  This  Table  refers  to  the 
superconvergence  points  given  in  Tables  la,  lb  for  the  class  of  harmonic  solutions. 
The  points  from  Tables  la,  lb  which  are  also  superconvergence  points  for  general 
solutions  for  this  pattern  are  indicated  by  the  word  yes. 

Table  4c.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  (1  <  p  <  7).  General  solutions:  Union-Jack  pattern.  This  Table  refers 
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to  the  superconvergence  points  given  in  Tables  la,  lb  for  the  class  of  harmonic 
solutions.  The  ue  no  superconvergence  points  for  general  solutions  for  this  pattern. 

Table  5a.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  (1  <  p  <  6).  General  solutions:  Criss-Cross  pattern.  This  Table  refers 
to  the  superconvergence  points  given  in  Tables  2a,  2b  for  the  class  of  harmonic 
solutions.  The  points  from  Tables  2a,  2b  which  are  also  superconvergence  points 
for  the  class  of  general  solutions  for  the  element  Ti  (shown  in  Fig.  5d)  are  indicated 
by  the  word  yes. 

Table  5b.  Superconvergence  points  for  periodic  meshes  of  triangular  elements  of 
degree  p  (1  <  p  <  6).  General  solutions:  Criss-Cross  pattern.  This  Table  refers 
to  the  superconvergence  points  given  in  Tables  3a,  3b  for  for  the  class  of  harmonic 
solutions.  There  are  no  superconvergence  points  for  the  class  of  general  solutions 
for  the  element  rj  (shown  in  Fig.  5d). 

Table  6a.  Superconvergence  for  meshes  of  triangular  elements.  Quadratic  el¬ 
ements,  Regular  pattern,  harmonic  solution  =  sin(frxi)sinh(wx3))  : 

Convergence  of  the  discrete  norms  for  various  selections  of  the  points.  Column  2: 
Convergence  of  the  discrete  norm  E'  which  includes  all  four  points  given  in  Table 
la  for  p  =  2.  Column  3:  Convergence  of  the  discrete  norm  E"  which  includes  only 
the  first  two  points  from  Table  la  (the  points  belong  to  the  element  bound' 
Column  4:  Convergence  of  the  discrete  norm  E'"  which  includes  only  the  third  and 
fourth  points  from  Table  la  (the  points  are  in  the  interior  of  the  element).  Note 
that  all  three  discrete  norms  E'"  are  superconvergent. 

Table  6b.  Superconvergence  for  meshes  of  triangular  elements.  Quadratic  ele¬ 
ments,  Regular  pattern,  general  solution  (u^(xi,X2)  =  sin(wxj)sin(wx3))  :  Con¬ 
vergence  of  the  discrete  norms  for  various  selections  of  the  points.  Column  2: 
Convergence  of  the  discrete  norm  E'  which  includes  all  four  points  given  in  Table 
la  for  p  =  2.  Column  3:  Convergence  of  the  discrete  norm  E"  which  includes  only 
the  first  two  points  from  Table  la  (the  points  belong  to  the  element  boundary). 
Column  4:  Convergence  of  the  discrete  norm  E'"  which  includes  only  the  third  and 
fourth  points  from  Table  la  (the  points  are  in  the  interior  of  the  element).  Note 
that  only  the  discrete  norm  E"  is  superconvergent. 

Table  6c.  Superconvergence  for  meshes  of  triangular  elements.  Cubic  elements. 
Regular  pattern,  harmonic  solution  («^(x,,X2)  =  sin(irxj)sinh(irx2))  :  Conver¬ 
gence  of  the  discrete  norms  for  various  selections  of  the  points.  Column  2:  Con¬ 
vergence  of  the  discrete  norm  E'  which  includes  all  six  points  given  in  Table  la  for 
p  =  3.  Column  3:  Convergence  of  the  discrete  norm  E"  which  includes  only  the 
third  point  from  Table  la  (the  mid-point  of  the  the  side  of  the  element  parallel  to 
the  x^-axis).  Column  4:  Convergence  of  the  discrete  norm  E'"  which  includes  all 
the  points  from  Table  la  except  the  third  one.  Note  that  all  three  discrete  norms 
E',  E",  E'"  are  superconvergent. 
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Table  6d.  Superconvetgence  for  meshes  of  triangular  elements.  Cubic  elements, 
Regular  pattern,  general  solution  =  sin(«'x J sin(irz3))  :  Convergence 

of  the  discrete  norms  for  various  selections  of  the  points.  Column  2:  Convergence 
of  the  discrete  norm  E'  which  includes  all  six  points  given  in  Table  la  for  p  =  3. 
Colunm  3:  Convergence  of  the  discrete  norm  E"  which  includes  only  the  third 
point  from  Table  la  (the  mid-point  of  the  the  side  of  the  element  parallel  to  the 
Xj-axis).  Column  4:  Convergence  of  the  discrete  norm  E'"  which  includes  aU  the 
points  from  Table  la  except  the  third  one.  Note  that  only  the  discrete  norm  E"  is 
supercon  vergent . 

Table  7.  Superconvergence  points  for  square-elements  of  the  serendipity  family. 
Superconvergence  points  for  the  x^ -derivative  for  the  class  of  harmonic  solutions  for 
elements  of  degree  p  (3  <  p  <  7).  The  vertices  of  the  square  are  at  (-1.,-!.),  (1.,-!.), 
(l.,l.),  (-1.,!.).  Note  that:  a.  For  p  =  3  we  have  three  lines  of  superconvergence  for 
the  class  of  harmonic  solutions;  the  thick  continuous  line  in  Fig.  5.  In  this  Table 
for  p  =  3  we  give  the  four  points  and  the  line  which  are  also  superconvergent  for 
general  solutions;  b.  For  p  >  4  there  are  no  superconvergence  points  for  general 
solutions. 

Table  8a.  Superconvergence  for  meshes  of  square  elements.  Cubic  elements,  gen¬ 
eral  solution  (u®(xi,X2)  =  sin(irxj)sin(irx2))  :  (3x3)  Gauss-Legendre  points. 
Values  of  the  discrete-norm  of  the  Xj-derivative  E  and  of  the  product  h~^E  for 
uniform  subdivisions  of  the  mesh-size  h.  Note  that  the  values  of  E  are  supercon¬ 
vergent  for  both  the  tensor-product  and  the  intermediate  space. 

Table  8b.  Superconvergence  for  meshes  of  square  elements.  Cubic  serendipity 
elements,  general  solution  (u^(x^,x,)  =  sin(7rxj)sin(irx2))  :  Effect  of  the  selection 
of  the  points.  Column  2:  Values  of  the  discrete  norm  of  the  x^-derivative  E  and 
of  the  product  h~^E  for  uniform  subdivisions  of  the  mesh-size  h  when  the  (3  x  3) 
Gauss-Legendre  points  are  used  in  the  computation  of  E.  Column  3:  Values  of 
the  discrete-norm  for  the  x^-derivative  E  and  of  the  product  h~*E  for  uniform 
subdivisions  of  the  mesh-size  h  when  the  points  given  in  Table  6a  are  employed  in 
the  computation  of  E.  Note  that  the  values  of  E  are  superconvergent  only  when 
points  from  Table  7  are  employed. 

Table  8c.  Super  convergence  for  meshes  of  square  elements.  Quartic  serendip¬ 
ity  elements,  points  from  Table  7  :  Effect  of  the  nature  of  solution.  Column  2: 
Values  of  the  discrete-norm  of  the  Xj-derivative  E  and  of  the  product  h^  ^E  for 
uniform  subdivisions  of  the  mesh-size  h  for  the  harmonic  solution  u^(Xj,X2)  = 
sin(7rxj)sinh(irx2).  Column  3:  Values  of  the  discrete-norm  of  the  Xj-derivative 
E  and  of  the  product  h~*E  for  uniform  subdivisions  of  the  mesh-size  h  for  the 
general  solution  tt^(x^,X2)  =  sin(irxj)sin(7rx2).  Note  that  the  values  of  E  are 
superconvergent  only  for  the  harmonic  solution. 
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Triangular  elements 


Superconvergence  points  for  the  class  of  harmonic  solutions 


The  vertices  of  the  triangle  are  at  (0. 


Regular,  Chevron,  Union-Jack  patterns 


Points 


.0118976218, 


tilililiTiIiTiTiIiTiTi 


0.3752181980, 

0.0858825088 

0.2821230756, 

0.0000000000 

0.5000000000, 

0.7178769244, 

■ilililililililiinilil 

■iIiIiItlilililililiM 

0.4715185566, 

E 

.9881023782, 

0.3241733208 

iiTiiiTiiiiiTiiiiiTin 


.1726731646, 
00 


iiimTiTiTiiiTiTliiTIi 


1.0000000000,  0.1195044522, 
0.0000000000  0.0224992688 

0.8638176260,  0.7053452590, 
0.6824066036  0.0071869378 

0.3095295948,  I 
0.2201118672 

0.6402277154, 
0.5196142358 

0.9063362308, 
0.7871241952 


liWiIiItlUtTiIilil 
[ULiIiIiIiIiIiIiIiM!! 
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Triangular  elements 


Superconvergence  points  for  the  class  of  harmonic  solutions 


The  vertices  of  the  triangle  are  at  (0 


Regular,  Chevron,  Union-Jack  patterns 


Points 


0.9871917636,  0.0797658369 
0.9556960448,  0.7771050728 


0.8901470515,  0.0000000000 


0.7965034496,  0.7034683949 
0.6451460051,  0.0000000000 
0.5873974006,  0.5156451535 
0.3658757490,  0.2978238991 
0.3548540007,  0.0000000000 
0.2851597869,  0.0102771287 


0.2807031259,  0.1264409904 


0.1098529492,  0.0000000000 


0.0637900495,  0.0000000000 
0.9837699874,  0.0278840787 
0.9671927841,  0.8195293891 
0.7443785569,  0.0000000000 
0.8413048661,  0.7026536445 
0.6689015074,  0.6017629175 
0.5000000000,  0.0000000000 
0.4770713108,  0.4203250119 
0.2949370204,  0.2365148149 
0.2556214434,  0.0000000000 
0.2236970856,  0.0897505908 
0.1620004419,  0.0056060808 
0.9362099504,  0.0000000000 
1.0000000000,  0.0000000000 
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Triangular  elements 


Superconvergence  points  for  the  class  of  harmonic  solutions 

The  vertices  of  the  triangle  are  at  (0.0,0.0) 

,  (1.0, 0.0),  (0.5,  0.5) 

Criss-Cross  pattern 

Points 

p=  1 

p  =  2 

■sa 

•e 

II 

■a. 

p  =  5 

1 

0.2113248654, 

0.0000000000 

0.1464466094, 

0.0000000000 

0.0761088279, 

0.0000000000 

0.0488705882, 

0.0000000000 

2 

0.7886751346, 

0.0000000000 

0.5000000000, 

0.0000000000 

0.3397707358, 

0.0000000000 

0.2355765633, 

0.0000000000 

3 

0.8535533906, 

0.0000000000 

0.6602292641, 

0.0000000000 

4 

0.9238911721, 

0.0000000000 

0.7644234366, 

0.0000000000 

5 

0.5000000000, 

0.4887686394 

0.9511294118, 

0.0000000000 

6 

■iiniTiTiiiiiiiTiiiTiil 

Table  2a 
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TrianguHar  elements 

Superconvergence  points  for  the  class  of  harmonic  solutions 

Criss-Cross  pattern 

Points 

p  =  6 

1 

0.4554029520,  0.4414423293 

2 

0.0344614509,  0.0000000000 

3 

0.5445970478,  0.4414423293 

4 

0.1716793301,  0.0000000000 

5 

0.3821587890,  0.0000000000 

6 

0.6178412109,  0.0000000000 

7 

0.8283206698,  0.0000000000 

8 

0.9655385491,  0.0000000000 

Table  2b 
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Triangular  elements 


Superconvergence  points  for  the  class  of  harmonic  solutions 


The  vertices  of  the  triangle  are  at  (O.OyO.O),  (1.0, 0.0),  (0.5,  0.5) 


Criss-Cross  pattern 


Points 

p=  1 

p  =  2 

i 

II 

p  =  4 

p  =  5 

■ 

1 

0.7886751346, 

0.5000000000 

0.8750000000, 

0.2834936534 

0.7621420044, 

0.3163540210 

0.8368698019, 

0.2268165908 

2 

0.8750000000, 

0.5000000000 

0.9145863153, 

0.3966634454 

0.7540881642, 

0.3397926069 

3 

0.8750000000, 

0.7165063710 

0.8970792676, 

0.5000000000 

0.6883014452, 

0.3781734456 

4 

0.5000000000, 

0.5000000000 

0.9145863153, 

0.6033365546 

0.7811733816, 

0.4582804499 

5 

0.7621420044, 

0.6836459789 

0.9334689552, 

0.5000000000 

6 

0.7811733816, 

0.5417195501 

7 

0.6883014452, 

0.6218265544 

i 

8 

1 

0.7540881642, 

0.6602073930 

9 

0.8368698019, 

0.7731834092 

10 

! 

ilM 

11 

0.7760422930, 

0.5000000000 

Table  3a 
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Triangular  elements 

Superconvergence  points  for  the  class  of  harmonic  solutions 

Criss-Cross  pattern 

Points 

p  =  6 

1 

0.7631964888,  0.2833600514 

2 

0.8508031659,  0.2069721253 

3 

0.8678258572,  0.1849536861 

4 

0.7086428573,  0.5000000000 

5 

0.8383266448,  0.3701490568 

6 

0.8284878782,  0.3953625943 

7 

0.8332694721,  0.500000000 

8 

0.8284878779,  0.6046374047 

9 

0.9467769341,  0.5000000000 

10 

0.8383266448,  0.6298509433 

11 

0.7631964887,  0.7166399487 

12 

0.8508031661,  0.7930278749 

13 

0.8678258572,  0.8150463139 

Table  3b 
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Triangular  elements 

Superconvergence  points  for  the  class  of  general  solutions 

wm 

Regular  pattern;  The  point  numbers  correspond  to  the  points  in  Tables  la,  lb 

Points 

a 

m 

m 

m 

B 

1 

yes 

yes 

no 

no 

no 

no 

no 

2 

yes 

no 

no 

no 

no 

no 

3 

no 

yes 

no 

yes 

no 

no 

4 

no 

no 

no 

no 

no 

no 

5 

no 

no 

no 

no 

no 

6 

no 

no 

no 

no 

no 

7 

no 

no 

no 

no 

8 

no 

no 

no 

yes 

9 

no 

no 

no 

no 

10  . 

no 

no 

no 

11 

no 

no 

12 

no 

13 

no 

14 

no 

15 

no 

Table  4a 
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Triangular  elements 


Superconvergence  points  for  the  class  of  general  solutions 

BB 

Chevron  pattern;  The  point  numbers  correspond  to  the  points  in  Tables  la,  lb 

Points 

B 

B 

B 

BB 

B 

wm 

IB 

1 

yes 

no 

no 

no 

no 

no 

no 

2 

no 

no 

no 

no 

no 

no 

3 

no 

yes 

no 

yes 

no 

no 

4 

no 

no 

no 

no 

no 

no 

5 

no 

no 

no 

no 

no 

6 

no 

no 

no 

no 

no 

7 

no 

no 

no 

no 

8 

no 

no 

no 

yes 

9 

no 

no 

no 

no 

10 

no 

no 

no 

11 

no 

no 

12 

no 

13 

no 

14 

no 

15 

no 

Table  4b 
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Triangular  elements 


Superconvergence  points  for  the  class  of  general  solutions 


The  vertices  of  the  triangle  are  at  (0. 


Union- Jack  pattern;  The  point  numbers  correspond  to  the  points  in  Tables  la,  lb 


Points 
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Triangular  elements 


Superconvergence  points  for  the  class  of  general  solutions 

Criss-Cross  pattern;  The  point  numbers  correspond  to  the  points  in  Tables  2a,  2b 

(For  the  element  tj  in  the  pattern) 

Points 

HI 

Bi 

wm 

B 

IB 

■B 

1 

yes 

yes 

no 

no 

no 

no 

2 

yes 

yes 

no 

no 

no 

3 

no 

no 

yes 

no 

4 

no 

no 

no 

no 

5 

no 

no 

no 

6 

no 

no 

7 

no 

8 

no 

Table  5a 
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Triangular  elements 

Superconvergence  points  for  the  class  of  general  solutions 

The  vertices  of  the  triangle 

Criss-Cross  pattern;  The  point  numbers  correspond  to  the  points  in 

(For  the  element  in  the  pattern) 

Tables  3a,  3b 

Points 

m 

la 

n 

m 

wm 

1 

no 

no 

no 

no 

no 

no 

2 

no 

no 

no 

no 

3 

no 

no 

no 

no 

4 

no 

no 

no 

no 

5 

no 

no 

no 

6 

no 

no 

7 

no 

no 

8 

no 

no 

9 

no 

no 

10 

no 

no 

11 

no 

no 

12 

no 

13 

no 

Table  5b 
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IViangular  elements 
Quadratic  elements,  Regular  pattern 
Harmonic  solution 


h 

E" 

E»i 

■ 

.06063740 

(3.88) 

.07027678 

(4.50) 

.04914226 

(3.15) 

IRI 

.00696365 

(3.57) 

.00806014 

(4.13) 

.00565854 

(2.90) 

.00084995 

(3.48) 

.00098347 

(4.03) 

.00069111 

(2.83) 

0.03125 

.00010559 

_ 

(3.46) 

.00012217 

(4.00) 

.00008587 

(3.81) 

Table  6a 
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Triangular  elements 
Quadratic  elements,  Regular  pattern 
General  solution 


h 

E' 

{h-^E‘) 

E" 

{h-^E") 

E«>  {h-^E"') 

.03445760 

(0.55) 

.03685704 

(2.36) 

.03187807 

(0.51) 

0.1250 

.00510295 

(0.33) 

.00494697 

(2.53) 

.00525431 

(0.34) 

0.0625 

.00087724 

(0.22) 

.00062918 

(2.58) 

.00106923 

(0.27) 

0.03125 

.00018505 

(0.19) 

.00007899 

(2.59) 

.00024949 

(0.26) 

Table  6b 
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IViangular  elements 
Cubic  elements,  Regular  pattern 
Harmonic  solution 


D 

£• 

E" 

[hr^E") 

E"'  {h-*E"') 

■1 

.00205242 

(0.53) 

.00184893 

(0.47) 

.00209075 

(0.54) 

.00012563 

(0.51) 

.00012009 

(0.49) 

.00012671 

(0.52) 

H 

.00000779 

(0.51) 

.00000750 

(0.49) 

.00000785 

(0.51) 

0.03125 

.00000049 

(0.51) 

.00000047 

(0.49) 

.00000049 

(0.51) 

Table  6c 
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Triangular  elements 
Cubic  elements.  Regular  pattern 
General  solution  u° 


h 

E' 

E" 

{h-*E*') 

E>« 

{h-^E") 

HI 

.01523855 

(0.98) 

.00172192 

(0.44) 

.01667522 

(1.07) 

.00211270 

(1.08) 

.00011551 

(0.47) 

.00231377 

(1.18) 

.00027112 

(1.11) 

.00000734 

(0.48) 

.00029697 

(1.22) 

0.03125 

.00003411 

(1.12) 

.00000046 

(0.48) 

.00003737 

(1.22) 

Table  6d 
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Square  elements 

Superconvergence  points  for  the  class  of  harmonic  solutions 

Points 

msm 

p  =  4 

p=  5 

mgm 

mggm 

1 

m 

-0.6587732224 

-0.7678487865 

-0.8542303405 

-0.7415311853 

-1.0000000000 

-1.0000000000 

-0.9584215542 

-1.0000000000 

2 

0.7745966692 

0.5687732224 

0.0000000000 

-0.5632551950 

0.5773502691 

-1.0000000000 

-1.0000000000 

-1.000000(K)00 

3 

-0.7745966692 

0.0000000000 

0.7678487865 

0.0000000000 

0.7415311853 

0.5773502691 

-0.5389584311 

-1.0000000000 

-0.9589718598 

-1.0000000000 

4 

-0.7745966692 

-0.6174062248 

-0.5773502692 

0.5632551950 

-0.8404899503 

-0.5773502691 

0.0000000000 

-0.7990568224 

-1.0000000000 

0.0000000000 

5 

0.7745966692 

0.6174062248 

0.5773502692 

0.8542303405 

-0.3896003525 

-0.5773502691 

0.0000000000 

-0.7990568224 

-0.9584215542 

0.0000000000 

6 

0.0000000000 

-0.5773502692 

-0.7861575208 

0.0000000000 

0.5389584311 

-0.1678536900 

0.0000000000 

0.0000000000 

7 

-0.5587732224 

0.5773502692 

-0.0538688187 

0.3896003525 

1.0000000000 

-0.1678536900 

0.0000000000 

8 

0.SS87732224 

-0.5494131405 

0.0538688187 

0.8404899503 

1.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

9 

0.0000000000 

0.7861575208 

-0.7415311853 

O.OOOOOOOODO 

0.0000000000 

1.0000000000 

10 

0.5494131405 

-0.8542303405 

0.0000000000 

0.0000000000 

-0.9584215542 

1.0000000000 

11 

-0.5773502692 

-0.5632551950 

0.7415311853 

0.1678536900 

1.0000000000 

1.0000000000 

12 

0.5773502692 

0.0000000000 

0.1678536900 

0.9589718598 

13 

-0.5773502692 

0.5632551950 

0.7990568224 

1.0000000000 

14 

0.5773502692 

0.8542303405 

0.7990568224 

0.9584215542 

15 

-0.7678487865 

1.0000000000 

16 

0.0000000000 

1.0000000000 

17 

0.7678487865 

1.0000000000 

Table  7 
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Square  elements 

(3  X  3)  Gauss-Legendre  points 

Cubic  tensor-product  space 

Cubic  intermediate  space 

D 

E  {h-*E  ) 

E  ih-*E  ) 

0.2500 

0.000691647  (0.177) 

0.00568511  (1.455) 

0.1250 

0.000043212  (0.176) 

0.00035669  (1.461) 

0.0625 

0.000002686  (0.176) 

0.00002229  (1.461) 

Table  8a 
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Square  elements 


Cubic  serendipity  space 

(3  X  3)  Gauss-Legendre  points 

Points  from  Table  7 

D 

E  (h-^E  ) 

E  (h-*E  ) 

0.2500 

0.05310222  (3.399) 

0.002463078  (0.631) 

0.1250 

0.00664258  (3.401) 

0.000156494  (0.641) 

0.0625 

0.00083252  (3.410) 

0.000009786  (0.641) 

Table  8b 


58 


Square  elements 
Quartic  serendipity  space  5^ 
Points  £rom  Table  7 


0.2500 

0.1250 


Harmonic  solution 

General  solution  u' 

1 

E  {h-*E  ) 

E  (h-*E  ) 

0.0625 


0.0011209  (1.148) 
0.0000350  (1.147) 
0.0000011  (1.148) 


0.0010429  (0.267) 
0.0000652  (0.267) 
0.0000041  (0.269) 
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Figure  la 


Figure  lb 


Figure  2b 


Figure  2c 


Figure  3b 


Figure  3c 


Figure  3d 


Figure  4a 


Figure  4b 


Figure  5 


Figure  6a 


Figure  6b 


Figure  6c 


Figure  6d 


Figure  7a 


Figure  8b 


Figure  8c 


Figure  8d 


Figure  8e 


Figure  8f 


Figure  8g 


Figure  8h 


Figure  9a 
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Figure  9c 


Figure  9e 


Figure  9f 


1 


Figure  9g 


The  Labonitoiy  for  Nunericol  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Technology  of  the  University  of  Mar^and,  under  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals; 

o  To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc,). 

Further  information  may  be  obtained  from  Professor  I.  Babu§lui,Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  CoUege 
Park,  Maryland  20742-2431. 


